TU-900 



(N 



O 

(N 



Making confining strings out of mesons 

Ryuichiro Kitano, Mitsutoshi Nakamura, and Naoto Yokoi 
Department of Physics, Tohoku University, Sendai 980-8578, Japan 



o 

< 

Abstract 
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■ The light mesons such as 7T, p, lu, jo, and ao are possible candidates of magnetic degrees 

of freedom, if a magnetic dual picture of QCD exists. We construct a linear sigma model 
to describe spontaneous breaking of the magnetic gauge group, in which there is a stable 
vortex configuration of vector and scalar mesons. We numerically examine whether such 
a string can be interpreted as the confining string. By using meson masses and couplings 
(ST as inputs, we calculate the tension of the string as well as the strength of the Coulomb 

force between static quarks. They are found to be consistent with those inferred from 
the quarkonium spectrum and the Regge trajectories of hadrons. By using the same 
Lagrangian, the critical temperature of the QCD phase transition is estimated, and a 
non-trivial flavor dependence is predicted. We also discuss a possible connection between 
^ ■ the Seiberg duality and the magnetic model we studied. 

S: 



1 Introduction 



The dual Meissner effect is a plausible explanation of the color confinement in QCD 0]. 
The condensation of the magnetic monopole, (m) 7^ 0, makes the QCD vacuum to be in the 
dual superconducting phase, where color fluxes sourced by quarks are squeezed into tubes, 
explaining the linear potential between quarks. 

In the QCD vacuum, there is another interesting phenomenon called the chiral symmetry 
breaking. It is believed that the quark-antiquark pair condenses in the vacuum, and the 
SU^Nf)^ x SU (Nf)n symmetry is spontaneously broken down to a diagonal SU (Nf)y group. 

Since two condensations, (m) 7^ and (qq) ^ 0, happen in the same dynamics, these 
two may be related. Indeed, lattice simulations of finite temperature QCD suggest that 
deconfinement and chiral symmetry restoration happen at similar temperatures HQ. This 
leads us to consider a simple unified picture: 

(qq) / Oh (mm) / 0, 

where "0" represents a non-abelian electric-magnetic duality. If the magnetic "monopoles" 
m and fh carry flavor quantum numbers, the condensations (m) = (m) 7^ describe Higgsing 
of the magnetic gauge group as well as chiral symmetry breaking. This phenomenon has been 
observed in supersymmetric gauge theories 

In this hypothesis, the "monopole" condensations give masses to magnetic gauge bosons 
and simultaneously provide massless pions as the Nambu-Goldstone bosons. There is in fact 
such a structure in the real hadron world. It has been known that the masses and interactions 
of the pions and the vector mesons such as the p and the uj mesons are well described 
by a spontaneously broken U(Nf) gauge theory Q. (See for an earlier discussion on 
the description of the vector mesons as gauge fields.) The vector mesons and the pions 
are respectively interpreted as the gauge fields and the uneaten Nambu-Goldstone bosons. 
Therefore, we are naturally lead to consider a possibility that the p and oj mesons are actually 
the magnetic gauge bosons of QCD. Along this line, it has been demonstrated recently that 
QCD regularized into an N = 1 supersymmetric theory has such a magnetic description 



via the Seiberg duality [11], where the magnetic Higgs fields m and m are the dual scalar 
quarks |l! 

The Higgs model of the U(Nf) gauge theory contains vector and scalar fields as well as 
strings as solitonic objects JT3-I16J. Since the string carries a magnetic flux in the magnetic 



picture, it can naturally be identified as the confining string via the electric-magnetic duality. 
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We examine in this paper whether such an identification works at the quantitative level. By 
using hadron masses and coupling constants as inputs, one can calculate the string tension and 
the Coulomb force between static quarks. We obtain values which are consistent with those 
inferred from the quarkonium spectrum and the Regge trajectories in the hadron spectrum. 

We write down a linear sigma model which includes the vector mesons, the pions and 
the scalar mesons in the next section. A vortex configuration in the model is constructed 
in Section |3|, and we compare the energy of the monopole-antimonopole system to the 
experimentally measured potential between a quark and an antiquark in Section ||. The 
critical temperature of the QCD transition is estimated in Section ||. The identification of the 
vector mesons as magnetic gauge bosons is motivated by recent discussions in supersymmetric 
gauge theories [12,17,18|]. (See also [19] for an earlier discussion.) We extend the discussion 
and propose a new interpretation in Section ||. 



2 Magnetic linear sigma model 

The magnetic picture of a confining gauge theory is supposed to be a Higgs model of some 
gauge theory. We apply this principle in QCD, and construct a model to describe Higgsing 
of the magnetic gauge group as well as chiral symmetry breaking. 



2.1 Lagrangian 



We propose the following Lagrangian to describe the magnetic picture of QCD. It is a 
U(Nf) gauge theory, and the Lagrangian possesses the U(Nj)l x U(Nf)n chiral symmetry. 
The vacuum expectation values (VEVs) of the Higgs fields, Hl and Hr, break the chiral 
symmetry down to the diagonal subgroup, U(Nf)v, providing massless Nambu-Goldstone 
bosons identified as pions and ij. The ry meson (or rf in the three-flavor language) can obtain 
a mass through a term which breaks axial U(l) symmetry explicitly such as det(HLH R ) 
although we ignore it in this paper. The VEVs of the Higgs fields give masses to U(Nf) 
gauge bosons. We identify these massive gauge bosons as the p and the u) mesonsQ. The 
Lagrangian is given by 

i i 

*~ 4 A"' r 4 r 

+^Tr [\D^H L \ 2 + \D^H R \ 2 ] 

-V(H L ,H R ). (1) 
*In the three-flavor language, one should include K* (892) and 0(1020) in the vector mesons. 
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The first and the second terms represent the kinetic terms of the U(l) and the SU (Nf) parts 
of the U(Nf) gauge bosons: ui^ and respectively. The Higgs fields Hl and H R are Nf x Nf 
matrices which transform as 



H L ->• giHLgjj 1 , H R ->• g H H R g R , 



(2) 



under the U(Nf)i J , the gauged U(Nf), and the U(Nf) R group elements, <7l, <7#, and c/r, 
respectively. The covariant derivatives are, therefore, given by 



D^Hl = d^H L + ig 2 H LP ^T a + igiQu^H^ 



(3) 



(4) 



Here, we normalized the SU(Nf) generators in the fundamental representation, T a , and the 
U(l) charge, Q, such that 



Tr (r a T b ^j =^S ab , 



and 



Q 



2N f 



(5) 



(6) 



The most general potential terms consistent with the symmetries are given by 
V(H L ,H R ) = fl 



A " A ! ( Tt{H l h[) + Tv(H R H R ) - 2 V ' 



+^ {Tr [{H{H L + H R H R f] - 4 (t±{H l h{) + Tr(H R H R )) } 
^-^(TT(H L Hl)-Tr(H R H R )) 2 

t TJ r _ 17 „ [ft \2 



+ 



A" 
+^Tr 



[H\H L - H R H R ) 



(7) 



where we assumed the parity invariance under Hl <h> H r . This potential stabilizes Hl and 
H R at 



(H L ) = (H R ) 



(8) 



At the vacuum, 4Nj degrees of freedom in Hl and ii^ break up to Nj massless Nambu- 
Goldstone bosons, Nj longitudinal modes of the gauge bosons, Nj massive scalar particles, 
and Nj massive pseudoscalar particles. The decay constant of the Nambu-Goldstone particles 
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is given by / w at tree level. The gauge group is completely broken and the unbroken global 
symmetry is vectorial U(Nf)y- 

The physical modes at the vacuum can be classified by the representations of U(Nf)y, 
the spin and the parity. The masses of the physical modes are given by 

singlet Nambu-Goldstone boson (77): = 0, (9) 

adjoint Nambu-Goldstone boson (ir): = 0, (10) 

singlet vector (w): m£ = g\fl, (11) 

adjoint vector (p): m 2 p = gjf 2 , (12) 

singlet scalar (/ ): m| = 2A /^, (13) 

adjoint scalar (do): m\ = 2\aJ%, (14) 

singlet pseudoscalar: m 2 PS = 2A ; f%, (15) 

adjoint pseudoscalar: vr?p A = 2X" f 2 , (16) 

at tree level. Terms with A' and A" are not very important in the following discussion^]. 
Hereafter, we take 

9i= 92 = 9, (17) 
as the p and uj mesons have similar masses. 



' The pseudoscalar particles are, in fact, CP even, and thus they are exotic states which are absent in the 
hadron spectrum. One should take large A' and A" to make the exotic states heavy so that the model can be 
a low-energy effective theory of QCD. We thank M. Harada, V.A. Miransky, and K. Yamawaki for discussion 
on this point. 
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2.2 Vector mesons and pions 



When we integrate out the massive scalar and pseudoscalar mesons, the model reduces to 
a non- linear sigma model of Ref. 0. The matching at tree level gives a = 1, where a is a 
parameter in the low-energy Lagrangian: 

C3 {1 ~f f * Tr[\d fl (U L U R )\ 2 ]. (18) 

The unitary matrices Ul and Ur are fields to describe the Nambu-Goldstone modes including 
the ones eaten by the gauge bosons. The transformation properties of Ul and Ur under the 
gauge and flavor groups are the same as Hl and Hr, respectively. From the low energy data, 



the preferred value of a is estimated to be a ~ 2 with an error of 15% [20|. Although there is 



a factor of two difference from the prediction, this discrepancy can be explained by including 



quantum corrections and/or higher dimensional operators. As discussed in Ref. [ 20 1 , the 
quantum correction makes the Lagrangian parameter a (A) approaches to unity when we take 
A to be large, such as a(A) ~ 1.33 ± 0.28 for A = 4:irf n ~ 1 GeV. Moreover, the quantum 
corrections from the scalar loops give positive contributions to the gauge boson masses, that 
further reduces the a(A) parameter. Therefore, one can think of the Lagrangian in Eq. (|l]) 
as the one defined at a high energy scale such as the mass scale of the scalar mesons. 

However, the large quantum corrections result in predictions which depend on the choice 
of input physical quantities when we work at tree level, although the differences should be 
canceled after including quantum corrections. In this case, one should choose a set of physical 
quantities which gives small enough coupling constants so that the use of the perturbative 
expansion is valid and the tree-level results are reliable. 

The Lagrangian has four parameters relevant for the discussion: g, f^, Ao, and A^. The 
Ao and A^ parameters can be obtained from the scalar masses as we discuss later. The gauge 
coupling constant g and the f v parameter can be estimated from two of physical quantities. 
The well-measured physical quantities which can be used as input parameters are [ pOfl : 

g p = (340 MeV) 2 , g pn7T = 6.0, F n = 92 MeV, m p = 770 MeV, (19) 

where g p and g p7Tn are the decay constant and the coupling to two pions of the p meson 
measured by p — > e + e~ and p — > ittt decays, respectively, and F w is the decay constant of the 
pion. The relations to the Lagrangian parameters at tree level are given by 

9p = gfl, Op™ = f , F n = fn, m p = gf n . (20) 
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Among them, the pair to give the smallest gauge coupling is g p and m p such as 

9 

9 P 



g=^ = 5.0, f n = 2£- = 150 MeV. (21) 



The value g = 5.0 means that the loop expansion parameter, g 2 N f / (Att) 2 , is of order 30% 
whereas other choices of input quantities give 90 — 210% for Nf = 2. Therefore, the choice 
above is unique to make a quantitative prediction. Indeed, the values in Eq. (21) are close 
to the ones evaluated at one- loop level. In Ref. f^(|, the parameters at a scale A ~ 1 GeV 
is obtained to be g(A) ~ 3.3 - 4.2, f w (A) ~ 130 - 150 MeV, and a(A) ~ 1.0 - 1.5, which 



reproduce all the physical quantities in Eq. (|19|), We use the values of g and /„. in Eq. (21) 
in the following discussion. However, we should bear in mind that there are theoretical 
uncertainties at the level of a factor of two in the results obtained at the classical level. 



2.3 Scalar mesons 

In the hadron spectrum, there are light scalar mesons, such as a, k, /o(980) and ao(980), 
which have not been understood as qq states in the quark model since they are anomalously 
light. We propose to identify them as the Higgs bosons in this linear sigma model. We do 



not consider heavier scalar mesons as candidates since otherwise the formulas in Eqs. (13) 
and (li~4]) indicate that the coupling constants are large and the perturbation theory would 
not be applicable. 

By taking the masses of /o(980) and ao(980) as input quantities^], i.e., 

m s = m A = 980 MeV, (22) 
Eqs. (FL3|) and (111) give the coupling constants Aq and Xa as 



V Ao = V Aa = 4.6, (23) 

at tree level, where f n in Eq. (|2l|) is used. We use these values of coupling constants for later 
calculations. 



3 Vortex strings 

Since the model has a spontaneously broken gauged £7(1) factor, there is a vortex string as a 
classical field configuration. The string carries a quantized magnetic flux. Below we construct 
a solution with a unit flux, which will be identified as the confining string. 

* Since a and k are quite broad resonances, we do not use their masses as inputs. 
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There have been similar approaches to the confinement in QCD. The Ginzburg-Landau 



models (the magnetic Higgs models) are constructed from phenomenological approaches [14, 
21, 22j or based on the QCD Lagrangian [23,p4| through the abelian projection |^5| , and 
the stable vortex configurations are identified as the confining string. In supersymmetric 
theories, there have been numbers of discussion on the vortex configurations 1(| 26-|31]. 
In particular, the non-abelian string jl5,16|, which we discuss shortly, has been extensively 
studied clS Bj candidate of the confining string. 

Our model combines Higgsing of the magnetic gauge group and chiral symmetry breaking. 
As discussed in the previous Section, the model parameters are fixed by physical quantities 
such as masses and couplings of hadrons. Therefore, the properties of the strings such as the 
string tension can be evaluated quantitatively. Below, we explicitly construct a classical field 
configuration of the vortex string. 



3.1 Non-abelian vortex solutions 

In this model, there are string configurations called the non-abelian vortices which carry the 
minimal magnetic flux. By defining the following gauge field, 

4=V^(Quy%+l^), (24) 

there is a vortex configuration made of, e.g., the i = j = 1 component rather than the overall 
U(l) gauge field uj^. Compared to the string solution made of w,, this non-abelian string 
carries only 1/Nf of the magnetic flux and thus it is stable. 

In constructing the vortex configurations, we follow the formalism and numerical methods 
of Ref. ]3^| , where the potential between a monopole and an anti-monopole is evaluated 
numerically in the abelian-Higgs model. Classical field configurations are constructed by 
numerically solving field equations while imposing the gauge field to behave as the Dirac 
monopoles |33| as approaching to their locations. 

We consider a non-abelian vortex solution, where the magnetic flux is sourced by a Dirac- 
monopole and a Dirac-antimonopole configurations of the gauge field with i = j = 
1, representing non-abelian monopole configurations. These monopole and anti-monopole 
are not present as physical states in the model of Eq. (||), and we introduce them as 
field configurations with an infinite energy, i.e., static quarks^ The object we construct 

§In U(N) gauge theories with Higgs fields in the adjoint representation, there are monopoles as solitonic 
objects which are identified as junctions of vortices (3^,Q rather than the endpoints. The monopoles we are 
considering should not be confused with such configurations. 



8 



here, therefore, corresponds to a bound state of heavy quarks such as the charmonium and 
the bottomonium. In order to describe light mesons, the light quarks should be present 
somewhere in the whole framework. We discuss a possible framework in Section ^. 

In the cylindrical coordinate, (p, <p, z), where the monopole and the antimonopole located 
on the z-axis at z = ±R/2, we denote (Ad)^ as the configuration to describe the monopole- 
antimonopole system. They are given by 

{A D yj = 0, (25) 
A 1 ^ = 0, except for i = j = 1, (26) 



and 

z-R/2 z + R/2 



All - Nflux 1 

A D = a D (p = 

V2g p 



(p. (27) 



\p 2 + (z- R/2) 2 ] 1 / 2 [p 2 + (z + R/2) 2 } 1 / 2 _ 
The number of the flux, N^ ux , is quantized as iVfl ux 6 Z by the Dirac quantization condi- 
tion Equivalently, the magnetic charge of the monopole is quantized as 

qm - (28) 

The gauge field is well-defined everywhere except for the interval —R/2 < z < R/2 on the 
z-axis. The Dirac quantization condition ensures that the interval is covered in a different 
gauge. For constructing a vortex configuration, the following ansatz are taken: 

A^ = A^, 4 = (A 2J )» + oJ t , (29) 

4 = 0, a 1 = a i {p,z)(p, (30) 

{Hh)ij = (H R )ij = (j>i{p, z)6ij , 4>i = 4%. (31) 
With the ansatz, the Lagrangian is reduced to 




(32) 
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and the field equations are obtained as 



(a 1 + a D S 



2N 



2N f 



(33) 



where we take the unit of 

A/w = i. 

For i ^ 1, a 1 = is the solution. 

The potential energy between the monopole and the anti-monopole is given by 

V(R) = 



(34) 



(35) 



9 2 R 
-!- / d 3 x 



9 ±2/ 1 , \ 1 

— ^(a + a D )a g ^ 



An 



JV 



/ 



'8JV* 



(36) 



The first term comes from the magnetic Coulomb potential, Vcouiomb = — Q'mag/^vrii. The 
second term is the contribution from the non-trivial field configurations, and gives the linear 
potential between a monopole and an antimonopole for a large R. The self-energies of the 
Dirac monopoles are subtracted, and thus this expression provides a finite quantity. 



For An = Xa, which is the case as in Eq. (23), the problem simplifies to the case of the 
abelian string. The field equations gives 



= 1, for % / 1, 

as solutions and the equations for <pi and a 1 becomes 



P J 2 



The potential energy is in this case given by 



V(R) 



g 2 R 



(37) 

(38) 
(39) 

(40) 



The Nf dependence disappears from the potential energy. 
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3.2 Numerical results 



We numerically solve Eqs. ( pq) and (f3^) by following the procedure explained in Ref. 32]. 
The partial differential equations are solved by using the Gauss-Seidel method. The obtained 
field configurations are used to evaluate the potential energy in Eq. (f40|). 

In the unit of Eq. (|35|), the potential energy V(R) times the gauge boson mass m p can be 
obtained as a function of m p R. In this normalization, we have a single parameter k defined 
by 

m s \/% \A~4 

k = —^ — = = • (41) 

V2m p g § 

This corresponds to the Ginzburg-Landau parameter of superconductors. The numerical 
results are shown in Fig. |], where the potential energies for iVfl ux = 1 are drawn with four 
choices of parameters, k = 0.1, 0.9, 1.7, and 2.5. We see a linear potential in a large R region. 
By fitting the slope of the linear regime, one can extract the string tension a in the unit of 
Eq. (H). We show in Fig. | the tension a as a function of k. These results are all consistent 
with Ref. [32 1, except that the unit of the flux is different due to the non-abelian feature 
of the vortex. For k = 1/V2, the field equations reduce to a set of first order differential 
equations whose solutions are known as the BPS state. In this case, the tension is simply 
given by a = tt, which we have confirmed with an accuracy of 0.1 — 0.2 percent. 



4 Comparison to QCD data 

Now we compare the numerical results with data from experimental measurements. We 
identify the non-abelian Dirac monopoles with the minimal magnetic charge, Nr ux = 1, as 
static quarks, since otherwise the string with iVfl ux = 1 is stable and such a stable string is 
absent in QCD. The potential between a quark and an antiquark with a distance R can be 
parametrized by the following form: 

V(R) = ~ + vR. (42) 

This potential, called the Cornell potential, well fits the quarkonium spectrum with param- 
eters: 

A ~ 0.25 - 0.5, y/a~430MeV. (43) 

A similar value of the string tension a is obtained from the Regge trajectories of the hadron 
spectrum. The lattice simulations also reproduce the shape of the potential with A ~ 0.25 — 
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0.4 [36-38J and yfa ' /m p ~ 0.50 — 0.55 [39|. In perturbative QCD, at tree level, the Coulomb 
part V ~ —A/R is obtained from the one-gluon exchange between quarks. At a higher loop 
level, the shape of the potential approaches to the form in Eq. (|42|) [40|. Computations at 



three-loop level have been performed recently in Refs. [41,42], and it is reported that the 



result is in good agreement with lattice simulations up to a distance scale R < 0.25 fm |41] 



See, for example, Ref. [43 for a review of the static QCD potential. 

The Cornell potential also well fits the numerically obtained potential in the previous 
section. We superimpose the fittings with the Cornell potential in Fig. |] as dashed lines. 



4.1 Coulomb potential 



In the electric picture, i.e., in QCD, the Coulomb part V ~ —A/R is obtained with 



A 



2N C 4vr' 



(44) 



where the strong gauge coupling g s depends on R through renormalization. 

By duality, in the magnetic picture, the Coulomb term is accounted by a magnetic 



Coulomb force between monopoles. By using the magnetic charge in Eq. (25) with 
the coefficient is given by 



b 



A = <L = 2tt 



47T 



g- 



(45) 
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This corresponds to the first term in Eq. (40). Using the value of g in Eq. (pi]), we obtain 

A = 0.25. (46) 

The value is consistent with Eq. (^). This is already an interesting non-trivial test of the 
hypothesis that the vector mesons are the magnetic gauge fields. 



Note here that the Coulomb term in Eq. (40) arises from a solution of the classical field 



equations with boundary conditions given by the Dirac monopoles. Although the vacuum 
is in a Higgs phase, the Coulomb force dominates when the distance R is small compared 
to the inverse of the gauge boson mass. In the world-sheet theory of the string, it has been 
known that the Coulomb force can be reproduced as the Luscher term which stems from 
the boundary conditions of the string world sheet ]44[| . Interestingly, the Luscher term gives 
A = n/12 ~ 0.26 which is pretty close to the above estimation. 

4.2 Linear potential 

As we have seen already, the linear potential is obtained as in Fig. ffl. The normalized string 
tension a is shown in Fig. |2]as a function of k. From Eqs. (|2l|), ( p?3| ) and (fil|), the k parameter 
is given by 

k = 0.90. (47) 
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With this value, we obtain from Fig. ^, 

a = 3.5. (48) 

to recover the mass dimension, we obtain 

yfa = 400 MeV. (49) 

This is close to \fa in Eq. (|43|). The prediction is not very sensitive to k. For example, 
K = 0.6 — 1.2 gives \fu = 360 — 420 MeV. Although we expect a large theoretical uncertainty 
from quantum corrections, it is interesting to note that the estimated string tension is in the 
right ballpark. The hypothesis that the p and uj mesons as magnetic gauge bosons and light 
scalar mesons as the Higgs bosons is found to be consistent with the experimental data. 

It is important to notice that there is no dependence on Nf in the field equations (|38|), 
(|39|) or in the expression of the QCD potential (f40|). It is essential to have this property that 
the string is non-abelian. The dimensionless quantity ^fa jm p is, in this case, predicted to be 
Nf independent, which is consistent with the results from the lattice QCD |39f . 



By using in Eq. (21 



5 QCD phase transition 



At a finite temperature, QCD phase transition takes place. The lattice simulations support 
that deconfinement and chiral symmetry restoration happen at similar temperatures. The 
chiral transition temperature has been computed in lattice simulations, and found to be 
T c ~ 150 — 160 MeV |3|,|4| for physical quark masses. 

A simple estimate of the transition temperature is possible in the magnetic model in 
Eq. ([!]). The deconfinement and the restoration of the chiral symmetry both correspond to 
the phase transition to the vacuum with Hl = Hr = 0, which is stabilized by thermal masses 
at a finite temperature. When we define the transition temperature T c to be the one at which 



the Higgs fields become non-tachyonic at the origin, the temperature is obtained to be [45 



T < = iw f f - (50) 



where the factor r/ is a dimensionless quantity given by 



, = 1 + gg + »»fe+"4 (51) 



m 2 s 3m| 



at the lowest level of perturbation. Each term in the r] parameter represents the contribution 
to the thermal masses of the Higgs fields from different particles. The first term, the unity, 
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is the contribution from the scalar mesons. One should add up all the particles which obtain 
masses from the VEVs of Hl and Hr. The estimation of rj is quite non-trivial since there 
are particles which we did not consider, such as nucleons, and also the summation should be 
weighted by the abundance in the thermal bath, which may be affected by their large thermal 
masses, i.e., there may be large higher order corrections. 

By putting f n in Eq. (|2l|), we obtain 

170 MeV x (|) , (N f = 2), 

V 3/ . ' (52) 



140 MeV x ( 1 ) , (N f = 3) 



.3, 

The value rj ~ 3 seems to give temperatures consistent with ones from lattice simulations. 
It is interesting that rj ~ 3 is obtained from Eq. (|5l| ) when we take mps around the cut-off 
scale, A ~ 1 GeV. 



The formula in Eq. (50) predicts that the transition temperature is inversely proportional 
to y/Nf. This is numerically consistent with the flavor dependence of T c studied in Ref. (3{| 
for two and three flavors in the chiral limit. There, T c is obtained to be 173 ± 8 MeV and 
154 ± 8 MeV for two and three flavors, respectively. A simulation with a larger number of 
Nf should be able to test this prediction. 

6 Non-supersymmetric duality from the Seiberg duality 

The assumption in the whole framework is the electric-magnetic duality between the SU (N c ) 
gauge theory with Nf massless quarks and U(Nf) gauge theory with bosonic Higgs fields. 
The replacement of N c in the gauge group with Nf is familiar in supersymmetric gauge 
theories. For example, the Seiberg duality in the N = 1 supersymmetric theories replaces 
SU(N C ) gauge group by SU(Nf — N c ) in the magnetic picture. We explain here a possible 
connection between the Lagrangian in Eq. ([!]) and the Seiberg duality, which is discussed in 



Ref. [12]. We extend the discussion of Ref. Jl2| regarding the vortex string and interpretations 
of constituent quarks. 

It is obvious that the non-supersymmetric QCD can be obtained from supersymmetric 
QCD's by adding masses to superpartners and send them to infinity. What is non-trivial is if 
a vacuum in the theory with small masses of superpartners is continuously connected to the 
non-supersymmetric theory when we send the masses to large values. Such a continuous path 
may or may not exist depending on the space of parameters defined by a supersymmetric 



theory to start with. Recently, it is found in Ref. [12] that there is an explicit model which 



15 





SU(N C ) 


SU(N f ) L 


SU(N f ) R 


U(1)b 


SU(N c )v 


U(1) B > 


U(1)r 


Q 


N c 


Nf 


1 


1 


1 





(N f - N c )/Nf 


Q 


K 


1 


N} 


-1 


1 





(N f - N c )/Nf 


Q' 


N c 


1 


1 





N~ c 


1 


1 


Q' 


N~ c 


1 


1 





N c 


-1 


1 



Table 1: Quantum numbers in the electric picture. 



reduces to QCD in a limit of parameters and has a vacuum with the same structure as the 
low energy QCD in a region of parameters where the Seiberg duality can be used. By hoping 
that the region extends to the QCD limit, one can study non-perturbative features of QCD, 
such as strings, at the classical level in the dual picture. 

The proposed mother theory is N = 1 supersymmetric QCD with iV c colors and Nf + iV c 
flavors. By giving supersymmetric masses to the extra N c flavors and soft supersymmetry 
breaking masses for gauginos and scalar quarks, one obtains non-supersymmetric QCD with 
N c colors and Nf flavors. The global symmetries and quantum numbers are listed in Table [l|, 
where SU(N C ) is the gauge group. The U(l) b> symmetry is absent in the actual QCD, and will 
be spontaneously broken in the vacuum we discuss later. In order to avoid the appearance of 
the unwanted Nambu-Goldstone mode associated with this breaking, we gauge U(1)b'- The 
SU(N c )v group is also an artificially enhanced symmetry, and thus we gauge it. Since the 
added gauge fields only interact with extra flavors, the limit of large mass parameters still 
gives the non-supersymmetric QCD we wanted. 

The magnetic picture of the mother theory is an SU(Nf) gauge theory with Nf+N c flavors 
and meson fields. The particle content and the quantum numbers are listed in Table It 
was found in Ref. [jl2| that there can be a stable vacuum outside the moduli space by the 
help of the soft supersymmetry breaking terms. The vacuum is at (q) = (q) ^ 0, where 
SU(Nf) x SU(Nj)l x SU(Nf)n is spontaneously broken down to a single vectorial SU(Nf)v 
symmetry, that is the isospin symmetry. The symmetry breaking provides massless pions 
and simultaneously gives masses to the SU(Nf) x U(1)b> gauge fields. Those massive gauge 
fields can be identified as the vector mesons, p and u. 

Although the deformation with massive iV c flavors provides us with a QCD-like vacuum, 



there are several unsatisfactory features as noted in Ref. [ 1 2 j . Here we discuss those issues 
and consider a possible interpretation. In the above discussion, it sounds somewhat strange 
that the U(1)b' gauge field is identified as the u; meson which is in the same nonet as the 
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1 


(27V/ - N c )/Nf 
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N f 
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-1 


(2N f - N c )/N f 



Table 2: Quantum numbers in the magnetic picture. 

p meson, whereas the U(1)b' seems to have a completely different origin from the SU(Nf) 
magnetic gauge group. Second, in the particle content in Table [2|, there are fields which 
have U(1)b charges ±1, i.e., "quarks." These degrees of freedom do not match the picture 
of confinement since they look like free quarks. Finally, there is a vortex string associated 
with the spontaneous breaking of U(1)b>, which we would like to identify as the QCD string. 
However, since the stability of the string is ensured by topology, it is stable even in the 
presence of the massless quarks. The real QCD string should be unstable since a pair creation 
of the quarks can break the string. 

A possible interpretation is emerged from the consideration of the origin of U(1)b' in the 
magnetic picture. As one can notice from the quantum numbers, U(1)b' in the electric and 
magnetic pictures look different. In particular, the gauged global symmetry in the electric 
picture is U(N C ) ~ (SU(N C ) x [/"(1))/Ztv c whereas one cannot find a U(N C ) gauge group in 
the magnetic picture. This leads us to consider a possibility that there is an additional U(l) 
factor as a part of the magnetic gauge group. The actual magnetic gauge group is U(Nf), 
and it is broken by a VEV of a field with the quantum number of Q' Nc q N f so that U(1)b> in 
the magnetic picture is an admixture of two J7(l)'s. Namely, the duality of the gauge group 
goes through an intermediate step: 

SU{N C ) x U(N C ) (electric) U(N f ) x U{N C ) (magnetic) 

-> SU(Nf) x SU(N C ) V x U(1) B > (magnetic). (53) 

Under this assumption, when we send the gauge coupling of f7(l)(c U(N C )) in the electric 
picture to be a large value, the gauge boson of the U(1)b> factor in the magnetic picture is 
mostly the one from the U(Nf) magnetic gauge group. The identification of the oj meson 
becomes reasonable since the origin is now the same as the p meson. 
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Since U (1) b' is spontaneously broken by (q) = (q) 7^ 0, there is a stable vortex string which 
can be explicitly constructed as a classical field configuration in the magnetic picture. The 
duality steps ( ^3|) imply that there is another string in the magnetic picture: one associated 
with U(Nf) and another with U(N C ). However, if we go back to the electric picture, there is 
only a single U(l) factor in U(N C ), which can only give a single kind of string. This sounds 
like a mismatch of two descriptions. 

We propose here that the U(Nf) string, made of q, q, p, and u, is in fact unstable since 
the "quarks" can attach to the endpoints, and thus that is the one which should be identified 
as the QCD string. The U(N C ) string is stable, but should decouple in the QCD limit. As 
mentioned already, there are "quarks" in the magnetic picture, q' , q' , Z and Z. They are 
natural candidates of the "quarks" which attach to the U(Nf) string. In turn, if they are the 
degrees of freedom at the string endpoints, a linear potential prevents them to be in the one- 
particle states. Therefore, the "quarks" disappear from the spectrum. This interpretation 
seems to give resolutions to all the unsatisfactory features raised before: the nature of u, free 
quarks, and the stable string. 

For this interpretation to be possible, q', q' , Z and Z should carry magnetic charges of 
U(Nf) in addition to the quantum numbers listed in Table H Since we assume the electric- 
magnetic duality between the SU(N C ) and the U(Nf) gauge groups, it is equivalent to say 
that q', q' , Z and Z should be colored under SU(N C ), i.e., Z and Z are the quarks (the 
non-abelian monopoles in the magnetic picture) and q' and q' are non-abelian dyons. It is 
interesting to notice that they indeed have N c degrees of freedom. 

In the SU(Nf) x U(1)b' magnetic gauge group, there is a U(l) factor which rotates a 
particular component of qi and qi, where I is the index of the SU(Nf) gauge group. The 
vortex string associated with such a U(l) factor is called the non-abelian string and the 
one with the minimal magnetic flux is stable. Therefore, the "quarks" should attach to 
this string. When we take q\ is the one which rotates under the U{1) factor and normalize 
the charge of it as unity, the charges of other charged fields are listed in the left column 
of Table H. By assuming that q and q have no magnetic charges, the Dirac-Schwinger- 



Zwanziger condition [46, 47 allows the magnetic charges listed in the right column of Table g 
as the minimal magnetic charges divided by (2ir/e) with e being the gauge coupling constant. 
Interestingly, they agree with the "color charge" of SU(N c )y up to a normalization, which 
may be indicating that a part of SU(N c )v in the magnetic picture descends from the 
electric gauge group, SU(N C ). For dynamical fields with both electric and magnetic quantum 
numbers, we loose the standard Lagrangian description of the model. However, since the 
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-1 
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-1 
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-1/JVc 


1 



Table 3: Electric and magnetic charges under a U(l) factor in SU(Nj) x U(1)b'- 

sector of q, q (and $) is all singlet under SU{N c )y and is decoupled from the colored sector, 
there can be a Lagrangian to describe it, and we assume that is the model in Eq. (|l|). 

It is amusing to see that many ingredients to describe the hadron world are present in this 
model, such as the vector mesons, the pions, the light scalar mesons, the QCD string, and 
the constituent quarks. This is somewhat surprising since the Seiberg duality is supposed 
to describe only massless degrees of freedom. The non-trivial success of the model may be 
indicating that the addition of N c massive quarks is a right direction to fully connect the 
electric and magnetic pictures of M = 1 super symmetric QCD. 
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